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Lewia Reaearch Center 
Cleveland, Ohio 


SUMMARY 


A multiple-grid convergence acceleration technique recently Introduced for 
application to the solution of the Euler equations by means of Lax-Wendroff 
algorithms Is extended to treat compressible viscous flow. 

Computational results are presented for the solution of the thin-layer 
version of the Navler-Stokes equations using the explicit MacComack algorithm, 
accelerated by a convective coarse-grld scheme. 

Extensions and generalizations are mentioned. 


INTRODUCTION 


N1 (1981) Introduced a multiple-grid convergence acceleration scheme for use 
with one-step Lax-Wendroff algorithms and applied It to the solution of the 
homoenthalplc Euler equations. Johnson (1982) extended this technique to the more 
general case of two-step Lax-Wendroff algorithms and Illustrated Its use. In 
conjunction with several members of this class, for the efficient solution of the 
full Euler equations. Of the two-step multiple-grid algorithms tested, all 
required substantially less work than Nl's original one-step method. The greatest 
efficiency %fas attained by accelerating the convergence of the MacComack (1969) 
algorithm. 

Although Initial applications have focused on the Euler equations, the 
multiple-grid scheme presently under discussion nay be employed with any 
hyperbolic system of conservation laws. Furthemore, It nay be used to good 
advantage In accelerating the convergence of time-asymptotic solutions to viscous 
flow problems. Since such viscous flow computations are quite time-consuming. 
Improvement In their efficiency Is necessary In order to facilitate their use as 
design tools. 

The present work discusses the extension of the multiple-grid convergence 
accelsratlon scheme for use In solving the Navler-Stokes equations. Two 
possibilities are described: the full coarse-grld scheme and the convective 

coarse-grld scheme. The full scheme Is a straightforward application of 
previously-described principles to the Navler-Stokes system of conservation laws. 
The convective scheme Is based on heuristic physlcsl arguments and results In a 


•lapliflad procadur*. 

Raaulta «r« prasantad for an application of tha convactlva coaraa-t^ld 
achaMt uaad with MacComack'a natkodt for tha aolutlon of tha thln-layar varalon 
of tha two-dlaanalonal Maviar-Stokaa aqiutlona. Tha choicaa to uaa tha thln-layar 
aquatlona and to work In two dlnanalona wara nada for convanlanca and do not iaply 
any raatrlctlona inharant In tha nultlpla-grld achaaa. Extanalona to thraa 
dinanalona and to tha full Navlar-Stokaa aquatlona ara conceptually clear and 
are raaarvad for future work. 


VISCOUS FLOW EQUATIONS 


At preaant, the efficient numerical almulatlon of compreaalble flowa by maana 
of model equatlona baaed on the aaaumed exlatence of a velocity potential function 
la quite routine. More general inviacid flow problema may be modelled uaing the 
Euler equatlona. Efficient computational procedurea for aolvlng theae aquatlona 
are maturing rapidly. A limited range of vlacoua flowa m^y be treated through a 
judicious combination of such inviacid models with correction procedures based on 
boundary layer theory. In the more general case, for example where the inter- 
action between inviacid and viscous flow regions is no longer weak, the usual 
boundary layer approximations cannot be used. Here, the approach providing the 
greatest versatility is the use of the full Mavler-Stokes equations, either in 
some sub-domain or in the entire flowfleld. 


Navier-Stokes Equations 

Following Peyret and Vlvland (1975), the nondimens Iona 1 Navier-Stokes 
equations may be written in conservation law form as 


q - -(f + g ) + Re"^(r + s ) 

^t X *y X y 

where : 
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H«r« p, u, V, p, « and E ar« raspactlvaly danalty, valocity coaponanta in tha x- 
and y- dlractlona* praaaura, aound apaad and total anargy par unit voluna. Tha 
total anargy par unit voluna nay ba axpraaaad aa 

B ■ p(a + ^(u^2)) 

whara tha apaclflc intamal anargy, a, la ralatad to tha praaaura and danalty by 
tha alnpla law of a calorlcally parfact gaa 


p - (Y-l)pe 


with Y denoting tha ratio of apaclflc beats. The coefficient of thermal 
conductivity, k, and the vlacoalty coafflclanta, X and y, are aaaumed to be 
functlona only of temperature. Furthermore, by Invoking Stokaa' aaaumptlon of 
aero bulk vlacoalty, X nay be expreaaed In tama of the dynamic vlacoalty y aa 



Re and Pr denote the Reynold a and Prandtl numbers, respectively. 

Although, for simplicity, the Navler-Stokea equations are presented here 
written In Cartesian coordinates, Vlvland (1974) has shown that their strong 
conservation law form may be maintained under an arbitrary tlne-<lependent 
transformation of coordinates. Explicit detail concerning the generalised 
coordinate version of these equations has been provided by Steger (1977) and 
need not be repeated here. 

Thin-Layer Equations 

The thin-layer approximation. In the words of Baldwin and Lomax (1978), "... 
evolves directly from a realistic assessment of %rhat Is really being computed In a 
typical high Reynolds number Navler-Stokea simulation." A highly stretched mesh 
Is used to resolve the large flow gradients normal to the vortlclty-generatlng 
surface. Consequently, because of limitations on cosiputer capacity, the diffusion 
terms Involving derivatives parallel to the surface are not resolved well enough 
to merit their computation. 

Similar viscous terms are also neglected In the classical boundary layer 
approximation. However, while the boundary layer approximation replaces the 
norma] momentum equation with the assumption that the normal pressure gradient Is 
sero across the viscous layer, all momentum equations are retained In the 
thin- layer approximation and no assumptions are made concerning the pressure. 
Consequently, the separation point Is not a singularity of the thin-layer model 
equations nor do the problssis associated with matching a boundary layer solution 
to an Invlscld outer flow occur when they are used. 
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In practlca, tha thln-layar aaauBptlon la laplanantad by ualng a body-flttad 
coordlnata ayatoa and naglacting tha vlacoua tana In tha coordinate dlractlon 
along the body. For Cartaaian coordlnataa, vlth z rapraaantlng the body- 
confomlng coordinate* the thin- layer version of Eqn.(l) is 


q. - “(f + g ) + Re”^ 8 
t z y y 


( 2 ) 


where : 
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The eztension of this assusptlon to generalized coordinates is straightforward 
and nay be found in Steger (1977). 


FINE-GRID SOLUTION ALGORITHM 


For convenience, Eqns. (1) and (2) nny be rewritten as 

q - -(F + G ) (3) 

z y 

where, for the full Navier-Stokes equations, 

F = f - Re ^r G = g - Re ^s 


while, for the thin-layer equations. 


F = f 


_l«v 

G = g - Re 8 


An approximate numerical solution to Eqn.(3) nay be computed by means of the 
two-step Lax-Wendroff algorithm Introduced by MacCormack (1969). Having chosen a 
grid representative of the resolution desired In the converged solution, the 
forward predictor - backward corrector version of MacCormack *s method may be 
written as 
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^ q(tH-At) - q(t) 




1.J 


First derivatives In the viscous terns ere bsckwsrd differenced in the predictor 
and forward differenced in the corrector. 


This approach to solving viscoi's flow problems is quite robust and has been 
in widespre^ and successful use for some time, both for the tlme-sccurate 
computation of unsteady flow and for the tlme-asyiq)totlc solution of steady flow 
problems. In the latter case, where accurate resolution of physical transients is 
not required, the numerical stability limitation inherent in this explicit method 
may severely restrict the speed of its convergence to the steady state. Providing 
a method to accelerate convergence in this case is the object of this report. 


COARSE-GRID ACCELERATION SCHEME 


Schemes of Lax-Wendroff type may be arrived at intuitively by using Taylor's 
theorem to write the approximation : 


At2 

q(t+At) - q(t) + Atq^ + -j- q^^ 


(4) 


Since we seek solutions to Eqn.(3), time derivatives may be expressed as space 
derivatives : 

’t - -<'x * 
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where A and B are the Jacobian matrices : 
A i 3F/oq B 3G/3q 

Substitution Into Eqn. (4) results In : 


q(t+-At) - q(t) - At(F +G ) 

X y 


At- 


A(F 4G ) 
X y 




(5) 


Second-order accurate spatial discretization of Eqn. (5) then yields a one-step 
Lax-Wendroff method. 
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For ttaapltt* v« aay make the folloviog finite-volume type epproxlmetlone : 


^1+1. J-1^ 


“ ^^1-1, j+1 ^^1-1, j ^1- 


8Ay ^®1-1,J+1 ”*■ “l,j+l ®l+l,j+l^ 

" ^®1-1, j-1 * ■*' ®l+l»j-l^ ^ 

If we define the "chec4(e" In q et cell centera auch that : 


j+1 ' “ Ax ^^1+l.J ^1+l.J+l^ ■ ^^1,J 'l.J+1^ J 

■ ^ ^®1.J+1 * ®1+1,J+1^ “ ^®1.J ^ ®l+l,j^ J 


It then follova that 


Consistent with the above approximations and definition, we may write the 
approximation : 

- At(P 4C ) 1 - 1 1 

^'*■2 ’^■*‘2 

This motivates the definitions : 

AFi l“^l 1^**1 1 


AGi 1"®1 l^**l 1 


We then make the approximations : 
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If we now define the "correction" to q at grid nodes : 
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Eqns. (6) and (9) constitute the one-step Lax-Wendroff method used as a basic 
Integration scheae by Ml (1981). He gives the following heuristic Interpretation 
to these equations : the first calculates the change In q occurring In a control 
volune during the Increment At while the second distributes the effects of the 
changes occurring In four nearest-neighbor control volunes to their coamon central 
nodal point where they are conblned to fom the correction to the vector of 
conservation variables* as Illustrated In Fig. 1. This Interpretation motivates 
the construction of the coarse-grld acceleration scheme. 
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Full Goaf -Grid SchM« 


Given the fine-grid corrections, which any be cos 4 >uted by sn> one- or 
two-step Lsx-Wendroff sch«M, ss shown in Job son (1982), we wish to use 
succv^ssively coerser grids to propsgste these corrections throughout the 
coaputationsl doasln, thus sccelerstlng convergence to the steady state while 
asinteining the accuracy deteralned by the fine-grid discretisation. Given s 
basic fine grid with the nuaber of points in each direction expressible ss 2^ -f 1 
for p s natural nuaber, let successively coarser grids be defined by successive 
deletion of every other point in each coordinate direction. The full coarse-grid 
acceleration scheae, as Illustrated in Fig. 2, then replaces the ccaputstion of 
coarse grid changes by Eqn.(6) with s restriction of the latest fine-grid 
correction. This restricted fine-grid correction is then distributed according to 
a coarse-grid version of Eqr. (9) to obtain a coarse-grid correction which is, in 
turn, prolonged to the fine grid to becoae the new fine-grid correction. One 
tine-cycle of the aultiple-grld scheae is coapoaed of an application of som 
Lax-Wendroff scheae on the fine grid followed by an application of the coarse-grid 
solution procedure to each successively coarser grid. The flow of Inforaation in 
this process is depicted in Fig. 3. 

In the basic integration scheme, a change at one grid point affects only its 
nearest neighbors while, in a k-level multiple grid scheae, the same change 
affects all points up to mesh spacings distant. Furthermore, since the 
change is always determined by Information from the fine grid and slaply 
propagated by the distribution formulae for coarser grids, fine grid accuracy 
is maintained. 

Convective Coarse-Grid Scheae 


The Bultiple-grid convergence acceleration procedure using the full 
coarse-grid scheme is essentially that described in Johnson (1982) for use with 
the Euler equations, with the Jacobian matrices cf the flux vectors replaced by 
their appropriate viscous counterparts. 

Consideration of the physical processes being modelled in a viscous flow 
computation leads to the form**' ^i.ion of an alternative coarse-grid scheme. 
Dissipative effects have a local character and their influence need not be taken 
into account in the construction of coarse-grid distribution foraulae. Rather, it 
is the convective terms, with their global character, which are the key eleaent in 
coarso-grid propagation. Hence, a coarse-grid scheae for viscous flew 
coaputations may be formulated on the basis of the Inviscld equations of motion. 
Such a convective coarse-grid scheae is inherently aore efficient than the full 
coarse-grid scheme because of the diminished computational effort associated with 
forming the Jacobian aatrices of the Euler flux vectors rather than those of the 
viscous flux vectors. An additional benefit is that the convective coarse-grid 
scheme leads to a multiple-grid convergence acceleration procedure which is 
Independent of the nature of the dissipative terms retained in the viscous model 
equations. That is to say that the coarse-grid scheme used previously for the 
Euler equations may be employed, without modification, to accelerate the 
convergence of viscous flow coaputations based on the Navier-Stokes equations, the 
thin-layer equations, or any other viscous model equations which contain the full 
inviscid Euler equations. 
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Th« corr«ctDMs of tho hourlotlc phyalcal romaonliig usod in th« foraulation 
of tho convoctivo cooroo-grld schoao to vorifiod by tho coaputotlonol rooulta 
proaontod balow. 


RESULTS 


Wa llluatrata the uae of nultiple-grid convergence acceleration in vlacoua 
flow compute tlona by meana of a mathematically and phyalcally almpla exe^>le. We 
uae MacCormack'a method aa a fine-grid aolutlon algorithm for the thin- layer 
veralon of the Navler-Stokea equatlona and accelerate Ita convergence to e ateacly 
aolutlon with the convective ccerae-grld achame. 

We conalder the aubaonic flow through a caacade of unataggered flat platea at 
xero angle of attack, aa llluatrated In Fig. 4. The ratio of exit atatlc pxcoaure 
to upatream total preaaure la 0.8430191, yielding flow Mach numbera In the 
vicinity of 0.5 for the teat caaea to be exhibited here. The Reynolda numbera, 
baaed on caacade gap and critical speed, span the approximate range from 8.4 X 10^ 
to 2.0 X 10^. Symmetry la Invoked to limit the alxe of the computational domain 
and the flow la aasumed to be Icmlnar. These assumptions are m^e for convenience 
In specifying the n\imber and location of fine-grid nodal points and do not Imply 
limitations on the generality of the method. 

The choice of boundary conditions Is also Indicated In Fig. 4. At the inlet, 
total temperature, total pressure and flow angle are fixed, while at the exit, 
only the static pressure Is specified. Along the plate surface the no-sllp 
condition Is applied and the temperature Is specified. Along the remaining 
boundary sections the tangency condition Is used. Uniform flow at the Isentroplc 
Mach minber implied by the ratio of exit static pressure to Inlet total pressure 
Is used as an Initial state. The values of the dependent variables on the domain 
boundax'les are updated only during the fine-grid computations. This decouples the 
coarse-grld acceleration scheme from the details of boundary condition 
implementation. 

As Illustrated in Fig. 5, three different fine grids are employed In this 
study. All have the same number of codal points and have their transverse grid 
lines located at the same positions. They differ In the positioning of their 
lateral grid lines. These are smoothly stretched away from the solid boundary In 
a geometric progression, starting from the Inltlnl spaclngs Indicated In Fig. 5. 
These fine grids each allow the construction of four successively coarser grids, 
as Indicated In Table I. The members of these grid families may then be used In 
combination to form multiple-grid sequences of length one through five. 

Computations have been performed for the combinations of Reynolds number and . 
fine grid configuration Indicated In Table II. Isomach contours for the converged 
solutions produced for each case are shown In Fig. 6. The contour levels 
displayed are not equally spaced and are the same for all five cases shown. 
Nevertheless, they provide a good qualitative Indication of the nature of the 
computed flowflelds. More quantitative Information Is contained In Fig. 7, where 
nomallxed u-veloclty profiles are Illustrated. Here, the u-velocity, normalised 
with Its value at the top boundary and same streaawlse station. Is plotted as a 
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function of rnlatlv* dlatnncn froa thn bottoa boundary. Curvoa for ovary aaeond 
straaMfiaa atation, atartlng with tha plata laading adga and andlng at tha outflow 
boundary, ara diaplayad. Thay ara ataggarad in proportion to tha apaclng of thalr 
raapactlva atraaawlaa atatlona. Tha data diaplayad In Flga. 6 and 7 ara for 
condltlona of optlnal work raductlon, as Indlcatad In Tabla II. Bowavar, In aach 
caaa tha aolution obtained la not a function of grid aaquanca length. Thla waa 
verified by axtenalve cootpu rational axpariaantatlon. 

Convarganca hlatorloa ara ahown In Fig. 8. For aach caaa, we dlaplay tha 
fine grid convergence hlatory and tha corraapondlng convarganca hlatory for that 
grid aaquenca length which produced tha beat work raductlon factor. Several of 
tha plota have been truncated to fit within the reaidual range diaplayad. 

For all five teat caaea and the five poaalble grid aequance langtha, tha 
computational work required to reduce a atandard error aaaaure to a apeclflad 
tolerance haa been recorded. In each caae, baaed on theae data, we may eatlnatc a 
work reduction factor and a corraapondlng optimal aultlple»grld aequence length. 
Here we define the work reduction factor to be that multiple by which the work 
required to produce a converged aolution ualng a alngle fine grid excaeda the work 
required to produce the aame result using that multiple-grid sequence length which 
minimizes the computational work. The results obtained are recorded in Tabla II. 
Work reduction factors ranging from 1.5 to 4.8 have been realized. We observe 
that while the work reduction factor and, possibly, the optimal grid sequence 
length decrease with Increasing grid stretching, they do not appear to decrease 
with increasing Reynolds number. 

Several Issues bearing on the multiple-grid work reduction factor should be 
mentioned at thla Juncture. When one considers convergence acceleration of tha 
turbulent full Navier-Stokes eqixatlcns, greater work raductlon than obtained hare 
will result by virtue of the Inclusion of the full viccous terms and turbulence 
modelling on the fine grid. The treatment of the Jacobian matrlcea used In the 
coarse-grld acceleration ncheme has a large Influence on the efficiency of the 
coarse-grld computations, and hence, on the work reduction factor. Substantial 
Improvement appears to be possible over the current treatment of these Jacoblans. 
In the present computations. Injection Is used as the restriction operetor and 
linear Interpolation Is chosen as the prolongation operator. These choices may 
not be optimal for use on highly stretched grids. Better choices could Increase 
both the optimal grid sequence length and the work reduction factor. Similar 
consequences might result from an alternative coarse-grid formation strategy. 

Given the encouraging results obtained to date, more comprehensive testing 
and more sophisticated applications of the viscous flow convergence acceleration 
Ideas presented here are planned. 


CONCLUSIONS 


TVo coarse-grld schemes for use In the multiple-grid convergence acceleration 
of viscous flow computations have been presented : the full coarse-grld scheme and 

the simpler and more efficient convective coarse-grld scheme. 
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Th« cuovACtiv* co«rs«-grld sehoM la of qulta gonarul api^licablllty and aaj 
ba uaad without ■odiflcation with any aodal aquatlooa in tha hiararchy ranging 
froa tha Bular aquationa to tha full Maviar-Stokaa aquationa. 

Coaputatlonal avidanca of tha utility of tha praaaut approach to convarganca 
accalaration haa baan provldad by uaing tha convactiva coaraa-grid achaaa* in 
conjunction with MacConack'a two-atap Lax-Wandroff aMthod, to aolva tha thln- 
layar varaion of tha Navlar-Stokaa aquationa for a ainpla nodal problan. 

Work raduction factora ranging fron 1.5 to 4.8 hava baan obtainad In Initial 
taating ovar a fairly broad range of Reynolds nuabara and grid atratchinga. 

Mora conprchanaive taating and wore aophlatlcatad appllcationa, including 
axtansion of tha tachniquea diacuased hara to throe dlmanalons, era plannad. 


UPDinCBS 


B«Idvin» B.S. and Lwmx, H.: Thin-Layar ApproxlBation and Alfabrale Modal for 

Saparatad T^irbulant Flawa. AIAA Papar 78-257, Jan. 1978. 

Johnaon, G.M.: M iltlpla-Grld Ar.calaratlon cf Laz-Uandroff Algorlttau. NASA 

1M-82843. Mar. 1982. 

MacCorvack, R.W.: Tha Effact of Vlacoaity in Byparvaloclty lapact Crataring. 

AIAA Papar ^ 9 . 354 , Apr. 1969. 

R.!i.: A Multlpla Grid Scheaa for Solving lha Eular Equations. AIAA Papar 

81-1025, Juna 1981. 

Payrat, R. and Vlvland, H.: Coaputatlon of Vlacouo Coapraaslbla Flows Baaed on 

tha Mavlar-Stokaa Equations. AGARD AC-212, Sep. 1975. 

Stager, J.L.: lapllclt Finite Difference Siaulatlon of Flow About Arbitrary 

Caoaatrlaa with Application to Airfoils. AIAA Papar 77-665> Juna 1977. 

Vlvland, H.: Foraas Conaarvatlvaa das Equations da la Dynamlqua das Gas. La 

Recherche Adrospatlala, No.l, Jan. -Feb. 1974, pp. 65-66. 


12 


TABLE I 


3RID DESCRIPTIONS 


Grid 

1 

2 

3 

4 

5 

Number of 
Points 

65 X 33 

33 X 17 

17 X 9 

9x5 

mmm 

mmm 


TABLE II. - ALGORITHM EFFICIENCY 


Test Case 

Reynolds 

Number 

Initial 

Transverse 

Spacing 

Optimal 

Sequence 

Length 

Work 

Reduction 

Factor 

a 


0.0125 

5 

4.8 

b 

8.4 X 10^ 

0.00625 

3 

2.7 

c 

3.4 X 10** 

0.00625 

2 or 3 

3.0 

d 

3.4 X lO** 

0.0025 

2 or 3 

1.5 







e 


2.0 X 105 


0.0025 


3 


1.8 
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Prolongation of Coarse-Grld 6q to Fine Grid 
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FIGURE 3. - Multiple-Grid Information Flow 
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FIGURE 4. - Physical Problem 
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Smallest Grid Spacing : Ax > 0.03125 

Ay - 0.0125 



Smallest Grid Spacing : Ax ■ 0.03125 
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Smallest Grid Spacing : Ax * 0.03125 

Ay ■ 0.0025 


FIGURE 5 


Fine Grid Configurations 
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FIGURE 6a. - Isomach Contours 

Re ■ 8. A X 10^ , Ay ■ 0.0125 



FIGURE 6b. 


Isomach Contours 

Re * 8.4 X 10^ , Ay ■ 0.00625 
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FIGURE 6c. - Isomach Cuntours 

Re - 3.4 X lO** , Ay - 0.00625 



FIGURE 6d. - Isomach Contours 

- 3.4 X lO** , 


Re 


Ay ■ 0.0025 


FIGURE 6e 



. - Isomach Contours 

Re - 2.0 X 10^ , Ay - 0.0025 
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FIGURE 8b. - Convergence 
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FIGURE 8c. - Convergence 
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FIGITIE 8d. - Convergence Histories 
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